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Abstract 

Equilibrium phase diagrams for neutron star matter at subnuclear densities are 
obtained at zero temperature. Spherical, rod-like and slab-like nuclei as well as 
spherical and rod-like nuclear bubbles are taken into account by using a compress- 
ible liquid-drop model. This model is designed to incorporate uncertainties in the 
nuclear surface tension and in the proton chemical potential in a gas of dripped 
neutrons. The resultant phase diagrams show that for typical values of these quan- 
tities, the phases with rod-like nuclei and with slab-like nuclei occur in the form 
of Coulomb lattice at densities below a point where the system becomes uniform. 
Thermal fluctuations leading to displacements of such nuclei from their equilibrium 
positions are considered through explicit evaluations of their elastic constants; these 
fluctuations can be effective at destroying the layered lattice of slab-like nuclei in 
the temperature region typical of matter in the neutron star crust. 
PACS : 26.60.-hc; 97.60.Jd 
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1 Introduction 



In the outer part of a neutron star, nuclei are considered to mainly determine the state 
of matter in equilibrium Except for a thin envelope of the star, due to Coulomb 
interactions, nuclei form a bcc lattice neutralized by a roughly uniform sea of electrons. 
With increasing density, weak interactions render nuclei neutron-rich via electron cap- 
tures. Then, at a density of about ~ 4 x 10^^ g cm~^, neutrons begin to drip out of these 
nuclei. The crystalline region of the star (not) including a gas of the dripped neutrons 
is usually referred to as an inner (outer) crust. In the deepest region of the inner crust, 
corresponding to densities just below the normal saturation density ps ~ 3 x 10^*^ g cm~^, 
not only are nuclei expected to have rod-like and slab-like shapes, but also the system is 
expected to turn inside out in such a way that the constituents of the original nuclei form 
a liquid containing rod-like and roughly spherical bubbles of the dripped neutrons. These 
transformations, as originally indicated by Ravenhall et al. and Hashimoto et al. 
stem from a delicate competition between the nuclear surface and Coulomb energies at 
small internuclear spacings. 

Recent calculations of the ground state of matter in the crust, performed by using 
specific nuclear models [||, |^, ||, indicate that at a density of order 10^^ g cm~^, rather 
small compared with ps, roughly spherical nuclei turn into elongated rod-like nuclei. Ac- 
cordingly, the lattice structure is changed from a regular bcc lattice to a two-dimensional 
triangular lattice. With further increasing density, these rod-like nuclei are transformed 
into slab-like nuclei, which are arranged in the form of a layered lattice. After that, a two- 
dimensional triangular lattice of rod-like bubbles and a bcc Coulomb lattice of roughly 
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spherical bubbles appear in turn. At last, at a density of about Ps/2, the system dissolves 
into uniform nuclear matter. Since slabs and rods look like "lasagna" and "spaghetti" , the 
phases with positional order of one- and two-dimension are often referred to as "pasta" 
phases. 

It is of interest to note the relevance of non-spherical nuclei to pulsar glitches and 
cooling of neutron stars. Pulsar glitches are usually considered to be manifestation of 
a sudden large-scale release of neutron vortices trapped by some pinning centers in the 
crust [0 . The force needed to pin vortices has yet to be clarified completely even for a bcc 
lattice of spherical nuclei due mainly to the uncertain properties of impurities and defects 
which may be involved in the vortex pinning p. Whether or not neutron vortex pinning 
occurs in the "pasta" phases possibly formed depends also on the unknown properties of 
impurities and defects in the lattice. The presence of non-spherical nuclei would accelerate 
the cooling of the corresponding region of neutron stars by opening semi-leptonic weak 
processes which are unlikely to occur for spherical nuclei [Q. This stems from the fact that 
in non-spherical nuclei, protons have a continuous spectrum at the Fermi surface in the 
elongated direction. However, such additional cooling would be rather small compared 
with the overall cooling of the stars, since the mass of the region where non-spherical 
nuclei are expected to occur is only a tiny fraction of the total mass of the star. 

At what densities the phases with non-spherical nuclei and with bubbles are energet- 
ically more favorable than the usual bcc phase and the phase of uniform nuclear matter 
depends on the properties of neutron-rich nuclei and of the surrounding pure neutron gas. 
The quantities that mainly describe such properties but are still uncertain are the nuclear 
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surface tension -Esurf and the energy required to add a proton to the pure neutron gas, i.e., 
the proton chemical potential iJ,^\ -Esurf controls the size of the nuclei and bubbles, and 
hence the sum of the electrostatic and interfacial energies. With increasing -Esurf and so 
this energy sum, the density at which the system becomes uniform is lowered. There 
is also a tendency that the lower the smaller p^. This is because —p^^ represents 
the degree to which the neutron gas favors the presence of protons in itself. This fea- 
ture is implied by the work of Arponen who studied the sensitivity of the properties 
of uniform nuclear matter at low proton fraction to the melting density p^. Although 
i^surf and p^^ have not been fully determined, previous authors (e.g., Refs. [|, |, |]) des- 
ignated the values of these quantities almost uniquely. In this paper, by generalizing a 
compressible liquid-drop model developed by Baym, Bethe and Pethick (hereafter 
denoted by BBP) in such a way as to incorporate uncertainties in -Esurf and p^\ we draw 
the equilibrium phase diagrams of zero-temperature neutron star matter at subnuclear 
densities. The resulting phase diagrams show that while the phases with rod-like bubbles 
and with spherical bubbles can exist only for unrealistically low -Esurf, the phases with 
rod-like nuclei and with slab-like nuclei survive almost independently of -Esurf and p^\ 
For these two phases, there are directions in which the system is translationally invari- 



ant. As noted by Pethick and Potekhin [|r^, this situation is geometrically analogous with 
a liquid crystal rather than with a rigid solid. Elastic properties of the nuclear rods and 
slabs are thus characterized by elastic constants used for the corresponding liquid- crystal 
phases, i.e., columnar phases and smectics A, respectively. They expressed such constants 
in terms of the interfacial and Coulomb energies. These energies, as well as the energies 
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of bulk nucleoli matter and of a electron gas, are well described by the zero-temperature 
approximation. This is because the temperature of the crust ^ 10^ K, as inferred from X- 
ray observations of rotation-powered pulsars, is small compared with typical nuclear and 
electronic excitation energies. However, such temperature may have consequence to the 
structure of the Coulomb lattices involved. For example, the outer boundary of the crust 
is acutely controlled by the temperature because of the density dependence of the melting 



temperature of the bcc lattice |12|. In the deepest region of the crust, thermal fluctua- 
tions may possibly spoil the ordering sustained by the layered and triangular lattices, as 
expected from the usual Landau- Peierls argument. We thus calculate thermally induced 
displacements of non-spherical nuclei by using the elastic constants given by Pethick and 
Potekhin |]1T[. The critical temperatures at which the ordered structure is destroyed by 



the thermal fluctuations are also estimated. 

In Section 2 we construct a compressible liquid-drop model for nuclei and bubbles, 
and we write down equilibrium conditions for zero-temperature matter in the crust. The 
resultant phase diagrams are displayed in Section 3. Properties of the phase transitions 
involved are discussed in Section 4. In Section 5 displacements of rod-like and slab-like 
nuclei are calculated at finite temperatures. Conclusions are given in Section 6. 

2 Model for neutron star matter at subnuclear den- 
sities 

In this section, we write the energy of zero-temperature matter in the crust and the 
conditions for its equilibrium. We consider five phases which contain spherical nuclei. 
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cylindrical nuclei, planar nuclei, cylindrical bubbles and spherical bubbles, respectively. 
Each phase is taken to be composed of a single species of nucleus or bubble at a given 
baryon density ri]^. Crystalline structures of these phases are taken into account in the 
Wigner-Seitz approximation, which leads to sufficiently accurate evaluations of the lattice 
energy In this approximation, a cell in the bcc lattice, including a spherical nucleus or 
bubble of radius tn, is replaced by a Wigner-Seitz cell defined as a sphere having radius 
Tc and the same center. A cylindrical nucleus or bubble having an infinitely long axis 
and a circular section of radius tn is taken to be contained in a cylindrical Wigner-Seitz 
cell having the same axis and a circular section of radius in place of a cell in the 
two-dimensional triangular lattice. For a planar nucleus with width 2rN, a cell in the 
one-dimensional layered lattice, having width 2rc, is identified with the corresponding 
Wigner-Seitz cell. The values of for these phases are chosen so that each Wigner-Seitz 
cell may have zero net charge and hence have the same volume as that of the original cell. 

2.1 Energy of matter 

For the purpose of constructing a formula for the energy of the system allowing for un- 
certainties in -Esurf and a compressible liquid-drop model for nuclei and bubbles is 
of practical use, which is characterized by the energy densities of uniform nuclear matter 
and neutron matter as well as by the interfacial and Coulomb energies. In this model, 
neutrons and protons are assumed to be distributed uniformly inside and outside the nu- 
clei or bubbles. A sea of electrons, which ensures the charge neutrality of the system, is 
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taken to be homogeneous. We thus write the total energy density -Etot as 



w^N + w^L + (1 - u)E^{n^) + Eciric) (nuclei) , 

W^N + + M-En(^n) + -E'e(?^e) (bubbles) 



where w-^ is the energy of the nuclear matter region in a cell as divided by the cell volume, 
wl is the lattice energy per unit volume as given for the finite-size nuclei or bubbles, is 
the number density of neutrons outside the nuclei or inside the bubbles, Ue is the number 
density of the electrons, E^ and E^ are the energy densities of the neutron matter and 
of the electron gas, respectively, and u is the volume fraction occupied by the nuclei or 
bubbles, 

d 



(nuclei) 



(bubbles), 

where d is the dimensionality defined as c? = 1 for slabs, d = 2 for cylinders and d = 3 for 
spheres, and n is the nucleon number density in the nuclear matter region. In Eq. (1), we 
ignore corrections due to nucleon pairing effects, shell and curvature effects in the nuclei 
and bubbles, and deformations and surface diffuseness of the nuclei and bubbles. 

In writing w^, we adopt a generalized version of the compressible liquid drop model 
developed by BBP ||10[. is then expressed as 



WN(n,x,nn,rN,rc,rf) 



un[{l — x)m^ + xmp]c^ + unW{k, x) 

+Wsurii'ri, X, nn, tn, u, d) + wc{n, x, r^, u, d) (nuclei) , 

(1 — u)n[{l — x)mn + xmp]c^ + (1 — u)nW{k^ x) 

+Wsuriin',x,na,rT^,u,d)+wcin,x,rT^,u,d) (bubbles) , 

(3) 



where i^p) is the neutron (proton) rest mass, W{k,x) is the energy per nucleon 
for uniform nuclear matter of nucleon Fermi wave number k = {3Tc'^n/2Y^^ and proton 



fraction x, as given by BBP [see Eq. (3.19) in Ref. ^^^], Wsurf is the nuclear surface energy 
per unit volume, and wc is the self Coulomb energy (per unit volume) of protons contained 
in a cell. 

W{k, x) includes four parameters characterizing the saturation properties of nearly 
symmetric nuclear matter, which are the binding energy wq, the saturation density 2A;Q/37r^, 
the incompressibility K and the coefficient s determining the symmetry energy. In the 
present work, the values for these parameters determined empirically by BBP in such 
a way as to reproduce nuclear masses and radii of /? stable nuclei are replaced by the 
following values that are generally accepted among recent literature: Wq = 16.0 MeV, 
ko = 1.36 fm^^, K = 230 MeV and s = 30.0 MeV. We remark in passing that such a 
replacement makes no significant difference in the phase diagrams for neutron star matter 
at subnuclear densities that will be exhibited in Section 3. 

The proton chemical potential fi^^ in pure neutron matter is contained in W{k,x). 
This quantity, dominating stability of the phases including a dripped neutron gas over 
the phase of uniform nuclear matter, has not been well determined from microscopic 
calculations based on nucleon-nucleon interactions. We thus set fx^^ (not including the 
rest mass) as 



/X. 



(0) 

p 



-C^nT , (4) 



where Ci is the free parameter being positive definite. As can be seen from Fig. 1, Eq. 
(4) approximately reproduces the overall density dependence of the results obtained from 



the Hartree-Fock theory with Skyrme interactions (see Ref. |]T3l and references therein) 
and from the lowest-order Brueckner theory with the Reid soft-core potential |jl^, |15|. 
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Hereafter, Cx will be set as Cx = 300,400,500,600 MeV fm^; the case of Ci = 300 MeV 
fm^ agrees well with the result of Siemens and Pandharipande [|I^ adopted by BBP, 
whereas the case of Ci = 400 MeV fm^ is consistent with the other results cited in Fig. 1 
for densities of interest to us n ^ 0.12 fm"'^. 

In constructing the surface energy per unit volume Wsmf, we take note of the expression 
given by BBP [see Eq. (4.30) in Ref. |jTO|] and the Hartree-Fock calculations by Ravenhall, 



Bennett and Pethick |]T6| (hereafter denoted by RBP) using a Skyrme interaction. For 
comparison, it is useful to consider the surface tension -Esurf that is independent of the 
shapes of nuclei and bubbles; this is related to Wsnri via 

'"N'U^surf /p-x 
-t^surf — } • lOJ 

ua 



Following the spirit of lida and Sato [17|, we express i?surf as 



E,^,, = C2tanhl^^jEfJ , (6) 

where C2 and C3 are the adjustable parameters that will be discussed just below, jj^^ = 
dE-^/dn^ — m^c^ is the neutron chemical potential in the neutron gas not including the 
rest mass, and 

BBP ^ ain-n,,f/^W{K,0)-W{k,x)] 
^"■•^ (367r) 1/3^0 ' 

with kn = (37r2n„/2)^/3 a = 21.0 MeV. If we set C2 = 1.0 and C3 = 3.5 MeV for 
Ci = 400 MeV fm^, as shown in Fig. 2, the surface tension calculated for matter in 
equilibrium (see Subsection 2.2 for details) agrees well with the RBP values for proton 
fractions of interest here x ^ 0.2. For Ci = 300,500,600 MeV fm^, consistency with the 
RBP values is also obtained from C2 = 1.0 and C3 = 3.5 MeV. Hereafter this value of C3 



will be fixed, whereas, for our present purpose of examining the dependence on nuclear 
models of the equilibrium phase diagrams for matter in the crust, we shall shake the 
values of C2 within a rather wide range covering C2 = 1.0 (see Section 3). 

We proceed to express the total electrostatic energy per unit volume wq+l = Wc + wi^. 
In the Wigner-Seitz approximation adopted here, it reads [|l] 

wc+L = ^^{exnr^Yufdiu) , (8) 

with 



Uu) 



d + 2 



d-2 \ 2 



u 



(9) 



The energy density of the pure neutron gas, whose density derivative together 
with -Esurf and fj.^^ , plays a role in determining the melting density pm llll , is set by using 

W{k,x) as 

^n(rin) = [W{K, 0) + mnc>„ . (10) 

Note that this energy density corresponds to the fitting formula determined by BBP [see 
Eq. (3.8) in Ref. |jlO|] so as to reproduce the data for ^ 1.7 fm~^ obtained by Siemens 
and Pandharipande in the lowest-order Brueckner theory with the Reid soft-core 



potential. These data hold good at least for ^ kg, as can be seen from comparison 
with the results from several works based on more elaborate many-body schemes and 
potential models [jl3|. We shall thus use the BBP formula without any modification. 

The energy density of the free electron gas, which is relativistic and degenerate, is 
finally expressed as 

3 

-Ee(ne) = -hckeUe , (U) 
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with the electron Fermi wave number ke = {Sn'^ricY^^ . Here, we ignore the electron mass 
and the electron Hartree-Fock energy. 

2.2 Equilibrium conditions 

As a next step towards the description of the equilibrium phase diagrams for crustal 
matter, we minimize the energy density i^tot with respect to five variables n, x, n^, 
and u under fixed baryon density rib, 

_ f un+{l-u)nn (nuclei) , , . 

1 {1 — u)n + un^ (bubbles) . 

and under charge neutrality, 

_ J xnu (nuclei) , , , 

^ 1 xn{l — u) (bubbles) . 

Such minimization results in four equilibrium conditions, which were given by BBP for 
spherical nuclei (see Section 2 in Ref. [|lOl). 

Optimization of E^^t with respect to r^r at fixed n, x, rin and u leads to a familiar 
relation for size-equilibrium, 

Wsurf = 2Wc+L • (14) 

This relation is thus independent of dimensionality d. Such independence was understood 
from the fact that the system, no matter what type of nucleus or bubble it may contain, 

is in a three-dimensional physical space. 

The condition for equilibrium of the dripped neutrons with the neutrons in the nuclear 

matter region, arising from minimization of E^ot with respect to n at fixed rib, nx, and 
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u, reads 



Here 



(15) 



nx,nn,rj^,u 



udn 

1 9(wn + ^l) 



1 — M (9n 



(nuclei) 



mnC^ (bubbles) 



(16) 



is the neutron chemical potential in the nuclear matter region, and 



rrinC + 



dUn 



rrinC + 



1 — u drin 

1 dwj^ 



(nuclei) 



(bubbles) 



(17) 



is the neutron chemical potential in the neutron gas modified by the surface energy. 

The condition that the nuclear matter region is stable against /3-decay is yielded by 
minimization of £^tot with respect to x at fixed n, rin, tn and u. It is expressed as 



where 



(18) 



(19) 



is the electron chemical potential, and 

9(wn + wl) 



ud{nx) 



nipC 



n{l-x),nn,r]si,u 



1 d{w^ + Wl) 



1 — u d{nx) 



— iripC 



(nuclei) 
(bubbles) 



(20) 



n(l-x),nn,rj<i,u 
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is the proton chemical potential in the nuclear matter region. 

Finally, minimizing E'tot with respect to tn at fixed un, x, Tc and (1 — w)nn for nuclei 
[(1 —u)n, X, Tc and un^ for bubbles] yields the condition for pressure equilibrium between 
the nuclear matter and the neutron gas, 



p(N) ^ p(G) 



(21) 



where 



p(N) 



du 

9(wn + wi,) 



nu,x,nii,rc 



du 



(nuclei) 
(bubbles) 



is the pressure of the nuclear matter region, and 



(22) 



(23) 



is the pressure of the neutron gas. 

At given riy, and subject to charge neutrality (13), we have calculated the total energy 
density i^tot for the five phases containing spherical nuclei, cylindrical nuclei, planar nuclei, 
spherical bubbles and cylindrical bubbles, respectively. In such calculations, we first 
express x, tn and Tc in terms of rib, n and nn by using conditions (12), (14) and (18). In 
this stage, we can rewrite A*n^\ jJ^^^ and P^^^ given by Eqs. (16), (17) and (22) as 



,,,,, , kdW(k,x) d (I d \ ^ 



(nuclei) , 



6 OK tn \i- — uj \n on J 



1 d \ 

+ Esuri + x[iie- {rrin - mp)c^] (bubbles) , 



(24) 



13 



^(^n, 0) + — + — (nuclei) , 

3 dfcn Tn 1 - M 



p(N) 



nkdW{k,x) d—1 dndEsmi 47rr 



2 

-E'surf + — " ^'^^^ + — ^(exn)^(l — m) (nuclei) , 



3 dk tn 9n d{d + 2) 

nA; (/c, x) d—1 dn u dE^uri 47rr 



+ ^i^^surf + - - ^^(exn)\l - u) (bubbles) 

ok tn tn 1 — 0*^ a(c( + 1) 

(26) 



Conditions (15) and (21) in turn give n and nn as function of Substitution of these 
values into Eq. (1) leads thus to the minimum values of i?tot at given nb and conformation. 
The energy density for /^-equilibrated uniform nuclear matter has also been evaluated at 
the same rib by using the energy densities nW{k, x) and Ee{ne) and by following a line of 
argument of BBP (see Section 8 in Ref. JlOll)- Here we ignore the presence of muons, the 
number density of which is, if any, far smaller than Ug for densities of interest here. We 
have finally found the phase giving the smallest energy density among the uniform and 
crystalhne phases; the resultant phase diagrams for the ground-state matter at subnuclear 
densities will be shown in the next section. 

3 Phase diagrams 

We may now draw the phase diagrams for neutron star matter in the ground state on the 
nb versus C2 plane for Ci = 300,400,500,600 MeV fm^.Q The results have been plotted 
in Fig. 3. We have thus confirmed the tendency, as mentioned in Section 1, that larger 

E'surf and lower fj,^^ help the uniform matter phase "erode" the phases with non-spherical 

""^Note that these diagrams are different from the usual ones written over the density versus temperature 
plane. 
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nuclei and with bubbles. 

As can be found from Fig. 3, the obtained phase diagrams basically reproduce the 
feature that with increasing density, the shape of the nuclear matter region changes like 
sphere —>■ cylinder —>■ slab — > cylindrical hole —>■ spherical hole —>■ uniform matter. This 
feature can be derived from simple discussions about the Coulomb and surface effects 
ignoring the bulk energy . We also find that the energy difference between two successive 
phases is generally of order 0.1-1 keV per nucleon (see Fig. 4); this result is consistent 
with that obtained by others in earlier investigations (e.g., Refs. 0, |^, |^). However, it is 
significant to notice that for most sets of Ci and C2 including the typical one (Ci = 400 
MeV fm^ and C2 = 1.0), the shape changes become simpler, i.e., sphere — > cylinder —>■ slab 
— > uniform matter. This behavior, due probably to bulk corrections to the Coulomb and 
surface effects, was not observed in the previous works (see, e.g., Refs. 0, ^ |[). Since the 
values of i?n, fJ'p^ and -Esurf as adopted by these works are consistent with those obtained 
here for Ci = 400 MeV fm^ and C2 = 1.0, we may conclude that the difference in the phase 
diagrams should originate mainly from that in the properties of uniform nuclear matter 
at low but finite proton fractions, i.e., x ~ 0.1. These properties depend on the adopted 
schemes for constructing interpolation between almost pure neutron matter and nearly 
symmetric nuclear matter in the absence of reliable microscopic calculations of the energy 
density of uniform nuclear matter at intermediate proton fractions. Better understanding 
of the equilibrium properties of neutron star matter at subnuclear densities will thus 
require future theoretical and/or experimental works covering the region of intermediate 
proton fractions. 
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Curvature effects in nuclei and bubbles, which were at least implicitly taken into 
account in Refs. 0, ^, ^ while have been ignored in the present analysis, may be another 
source of the difference in the phase diagrams from such earlier works. This is due to the 
conformation dependence of these effects: t^curv ( sphere) > Wcurv( cylinder) > Wcurv(slab) 
= > lUcurv (cylindrical hole) > Wcurv( spherical hole), where Wcnw is the curvature energy 
per unit volume. The degree of the resultant stability of the higher- density crystalline 
phase over the lower-density one depends on the values of Wcun, which are also uncertain 
at X ~ 0.1. If we determine these values from the curvature thermodynamic potential 
obtained by Kolehmainen et al. |[18[ in the Thomas- Fermi theory with Skyrme interactions, 
the energy gain of the higher-density crystalline phase with respect to the lower-density 
one amounts to a few keV per nucleon. For our typical nuclear model (Ci = 400 MeV 
fm^ and C2 = 1.0), as can be seen from Fig. 4, such curvature effects do not lead to 
appearance of the density region where the phases with bubbles are stable, but do lower 
the transition density from sphere to cylinder and that from cylinder to slab by 10% or 
so. In Fig. 3, correspondingly, the region surrounded by the bcc phase and the uniform 
phase is expected to enlarge. 

4 Properties of phase transitions 

The structural transitions examined in the previous section are first-order transitions. 
This is because the above-mentioned sequence of geometrical structure is mainly deter- 
mined by the competition between the surface and Coulomb energies. In order to confirm 
this property, we have plotted in Fig. 5 the sizes of the nucleus or bubble and of the 
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Wigner-Seitz cell, tn and rc, evaluated as functions of rih for Ci = 400 MeV fm^ and 
C2 = 1.0. It can be clearly seen from this figure that and Tc jump at the transition 
points. It is of interest to note that for C2 = 0.01, in which all the types of crystalline 
phases occur, the cell size shows a dependence on the dimensionality d such that it 
is largest for d = 3 (spheres and spherical holes) and smallest for d = 1 (slabs). This 
behavior was also obtained from the previous liquid-drop- model calculations 0]. 

The neutron density profile is useful in probing how the system dissolves into uniform 
matter. In Fig. 6 we have thus plotted the neutron densities in the nuclear matter region 
and the neutron gas region for the crystalline phases as well as the one for the uniform 
phase, evaluated as functions of nb for Ci = 300,400,500,600 MeV fm^ and C2 = 1.0. 
For these parameter sets, as shown in Fig. 3, there is a melting transition from the phase 
with slab-like nuclei to the uniform phase. We can observe from Fig. 6 that as the 
density approaches the melting point, the neutron distribution becomes more and more 
smooth. This behavior, as we have confirmed for various values of C2, is consistent with 
the results from the Thomas- Fermi calculations |T^. We also find that at the melting 
point, the neutron density profile is discontinuous between the crystalline and uniform 
phases, indicating that the transition is of first order. However, this discontinuity should 
be fairly small as compared with the case of matter created in stellar collapse. In this 
case, due to larger proton fraction x yielded by degenerate neutrinos, such a discontinuity 
should become so large that neighbouring nuclei (or bubbles) appear to touch and fuse 
with each other at the melting point. Detailed estimates made for this material will be 
reported elsewhere. 
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5 Thermal fluctuations 



Let us now estimate thermally induced displacements of rod-like and slab-like nuclei from 
their equilibrium positions; the presence of these nuclei is energetically favored for typical 
values of Ci and C2, as shown in Section 3. For such estimates, we first write the elastic 



constants for the phases containing these nuclei by following Pethick and Potekhin |[TT 



Since the phase with slab-like (rod-like) nuclei is structurally similar to a smectic A 
(columnar phase) liquid crystal, they derived the expressions for the elastic constants by 
comparing the energy increase due to deformation, obtained from the liquid-drop model 
for the nuclei in the incompressible limit, with that used for the corresponding liquid 
crystals. For the layered phase composed of slab-like nuclei, the energy density due to 
displacement v of layers in their normal direction (taken to be parallel to the 2;-axis) can 
be written as 

(27) 



2 



9V 1 ,n 



2 



+ ^(Vlt;)^ 



Here, the elastic constants B and Ki are expressed as [11 



B 



2 
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u;c+l(1 + 2u- 2uy, 



c 



(28) 
(29) 



For the two-dimensional triangular lattice of rod-like nuclei, the energy density due to a 
two-dimensional displacement vector v = {vx,Vy) of cylinders running along the z-axis is 
given by pO 



B ( dvr. dv.,, 



2 \dx dy 



~2 



' dvx _ dvy 
dx dy 



dy dx 



Here, the elastic constants B, B', B", C and read |11 



B = ^wc+L, (31) 

B' = -^wc+L, (32) 

B" = Iwc+L, (33) 

C ~ 102-^("-°-=^)wc+L , (34) 

~ 0.0655wc+L (for 0.15 <u< 0.55) . (35) 

In Eqs. (28), (29) and (31)-(35), wq+l, as given by Eq. (8), satisfies the size-equihbrium 
condition (14). 

We may then calculate the mean-square displacement at finite temperature T, defined 
as where (■ ■ ■) is the average over the probability distribution proportional to 

exp{— Jy dV F / k-QT) with the lattice volume V and the Boltzmann constant k^-^ Let us 
assume that the length scale L of the lattice is far larger than r^. For the one-dimensional 
layered lattice, the mean-square displacement can then be evaluated within the harmonic 
approximation allowing for the terms up to 0{v^) in Eq. (27) as PT 



, (36) 



^Hcre the fluctuations are treated thermodynamically. At T = 0, however, the fluctuations behave 
quantum-mechanically; the resultant mean-square displacements are of order Tlc2 / (r'^wc+h) , where C2 is 
the velocity of second sound associated with the varying inter- nuclear distance roughly estimated 
as C2 ^ \/wc+iJpm- For typical values rc ^ 10 fm, wc+l ~ 10~^ MeV fm~^ and p,„ ^ 10"^'' g cm~^, 
comparison of the quantum fluctuations thus estimated with the thermal fluctuations to be obtained 
below gives rise to the quantum-classical crossover temperature Tq of order lO'^ K for planar nuclei and 
of order 10^ K for cylindrical nuclei. Consequently, the forthcoming arguments hold good for T ^ Tq but 
underestimate the displacement for T ^ Tq. 
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where a = 2rc is the layer spacing. For the two-dimensional triangular lattice, by retaining 
the terms up to O(f^) and 0{Vy) in Eq. (30), we obtain 

^ ~ iB + 2C)V^ ' ^^^^ 

where A = \J'2,K^/ {B + 2C), and a = {2tt / ^/sY^'^Vc is the lattice constant for hexago- 
nal cells.0 Expressions (36) and (37) are analogous to those appearing in the context 
of Landau-Peierls instabilities, which destroy a one-dimensional ordering in an infinite 
three-dimensional system. 



In Fig. 7 (8) we have plotted the ratio y (|t;|^)/(a/2 — r^) of the root-mean-square 
displacement of the slab-like (rod-like) nucleus, calculated from Eq. (36) [(37)] for /cbT = 
0.1 MeV, to the shortest distance between the surface of the nucleus in its equilibrium 
position and the boundary of the cell containing it. In these calculations we have set 
L = 10 km, which is the typical neutron star radius. The results for the layered phase, 
depending on L logarithmically, are essentially unchanged by the choice of L; when we 
set L = 1 yum, for example, the values shown in Fig. 7 are only multiplied by a factor 
of ~ 0.7. It is obvious from comparison between Figs. 7 and 8 that for the same Ci 
and C2 the relative displacement of the slab-like nucleus is appreciably large compared 
with that of the rod-like nucleus, a feature yielded by the logarithmic factor appearing 
in Eq. (36). We also find that the relative displacement for both lattices decreases with 
increasing surface tension (or C2). This is because the elastic constants, acting to reduce 

the displacements given by Eqs. (36) and (37), are proportional to the equilibrium value 

■^Strictly speaking, Eq. (37) is valid when the length of the nuclei is large compared with the linear 
dimension of the lattice in the xy plane. This situation is energetically preferred over the opposite case, 
in which there is a larger total sectional area of the nuclei at the boundary of the lattice. 
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of wc+L satisfying condition (14). 

It is noteworthy that for the typical parameter set (Ci = 400 MeV fm^ and C2 = 1-0), 
the relative displacement for the layered phase at k-^T = 0.1 MeV, as plotted in Fig. 7, 
takes on values of up to about unity. The planar nucleus deformed to such an extent 
touches the boundary of the neighbouring cell as defined in the case in which the lattice 
is in its equilibrium. This suggests the possibility that the fluctuational displacements 
destroy the ordered configuration in a length scale of down to ~ Vc even at temperatures 
typical of matter in the neutron star crust .0 In theoretically predicting the interiors 
of observed neutron stars having various temperature profiles, therefore, it is useful to 
estimate as a function of rib the critical temperature Tc at which the relative displacements 
become unity. The results obtained for the layered lattice and the triangular lattice have 
been plotted in Figs. 9 and 10, respectively.^ The difference in between these two 
lattices, as can be observed from Figs. 9 and 10, suggests that if formation of these two 
lattices can occur dynamically in the star, the layered phase is formed later than the 
triangular phase during the star' s cooling. 

6 Conclusion 

We have examined the dependence on the surface tension -Esurf and on the proton chemical 

potential fj,^^ in pure neutron matter, of the density region in which the presence of non- 

^ Reconnection of the planar nuclei that may be involved in such destruction leads inevitably to 
departures from thermodynamic equilibrium, which are beyond the scope of this paper. 

^ Note that the values of Tc are meaningless when they are smaller than the crossover temperature Tq 
from the classical to the quantum fluctuations; in this case, quantum effects keep the relative displacement 
large compared with unity for T < Tq. 
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spherical nuclei and of bubbles is energetically favored at T = 0. We have found that as 
Esuri decreases or fi^^ increases, such a density region becomes larger. For the values of 
i^surf and fi^^ as adopted in recent literature, our results show that in the ground state, the 
phases with rod-like nuclei and with slab-like nuclei lie between the bcc lattice phase and 
the uniform nuclear matter phase. The fluctuational displacements of such non-spherical 
nuclei from their equilibrium positions have been estimated at finite temperature. It has 
been suggested that at temperatures typical of matter in the neutron star crust, such 
fluctuations may melt the layered lattice of slab-like nuclei. 

Even if rod-like and slab-like nuclei are thermodynamically stable, whether or not they 
are actually present in the star depends on occurrence of the dynamical processes leading 
to their formation. These processes are thought to involve instabilities against nuclear 
deformation and against proton clustering in uniform nuclear matter [0]. If the occurrence 
of such processes were confirmed, it would become still more significant to consider effects 
of the presence of non-spherical nuclei on the structure and evolution of neutron stars. 
In the context of pulsar glitches, the question of what kind of impurities and defects are 
formed in the lattices of non-spherical nuclei would be essential to understanding of the 
mechanism for neutron vortex pinning. For neutron star cooling, it would be interesting 
to consider the direct URCA processes in non-spherical nuclei; as suggested by Lorenz 
et al. these processes might be allowed by the proton continuous spectrum in the 
elongated direction of the nuclei, in contrast to the case of roughly spherical nuclei, and 
by the band structure of neutrons moving in the periodic potential created by the nuclei. 
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Fig. 1. The proton chemical potential in pure neutron matter as a function of neutron 
density. The thick lines are the results calculated from Eq. (4) for Ci = 300, 400, 600 MeV 
fm^. The thin broken lines as marked by the Skyrme interactions (FPS21, 1', EPS and 
SkM) are the results summarized by Pethick, Ravenhall and Lorenz |T^, and the thin 
solid line is the result of Sjoberg ||l5l. The crosses denote the values obtained by Siemens 
and Pandharipande . 



Eig. 2. The surface energy per unit area (the surface tension) as a function of x, the 
proton fraction in the nuclear matter region. The thick broken curves are the present 
results obtained from Eq. (6) for Ci = 400 MeV fm^ and Ca = 0.01,0.1, 1.0,2.5,5.0, the 



solid curve is the RBP result from their Hartree-Eock calculations [0, and the dotted 



curve is the BBP result ITO 



Eig. 3. Zero-temperature phase diagrams on the rib versus C2 plane, evaluated for 
Ci = 300, 400, 500, 600 MeV fm^. 

Eig. 4. Energy per nucleon relative to that of uniform matter calculated for Ci = 400 
MeV fm^ and C2 = 1.0 as a function of baryon density ny,. The symbols SP, C, S, CH 
and SPH stand for sphere, cylinder, slab, cylindrical hole and spherical hole, respectively. 

Eig. 5. Size of a nucleus or bubble, tn, and of a Wigner-Seitz cell, r^, calculated for 
Ci = 400 MeV fm^ and C2 = 0.01, 1.0, 2.5 as a function of baryon density Uh- The symbols 
SP, C, S, CH and SPH stand for sphere, cylinder, slab, cylindrical hole and spherical hole, 
respectively. 

Eig. 6. The neutron densities obtained for Ci = 300,400,500,600 MeV fm^ and 
C2 = 1.0 as functions of baryon density n\y. The lines classified by [n{l — x)]n and 
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represent the neutron densities in the nuclear matter region and in the neutron gas region 
for the phases with nuclei, respectively. The lines classified by [n{l — x)]uni denote the 
neutron densities in uniform nuclear matter. 

Fig. 7. The root-mean-square displacement of a planar nucleus at /cbT = 0.1 MeV, 
divided by the shortest distance between the surface of the nucleus in its equilibrium 
position and the boundary of the cell containing it. The curves are obtained for various 
sets of Ci and C2 as a function of baryon density 1%^. The thick curves lying between 
the two vertical lines are the results in the density region in which the phase with planar 
nuclei is energetically stable. 

Fig. 8. The root-mean-square displacement of a cylindrical nucleus at /cbT = 0.1 MeV, 
divided by the distance between the surface of the nucleus in its equilibrium position and 
the boundary of the cell containing it. The curves are obtained for various sets of Ci and 
C2 as a function of baryon density rih. The thick curves lying between the two vertical 
lines are the results in the density region in which the phase with cylindrical nuclei is 
energetically stable. 

Fig. 9. The critical temperature Tc for the phase with planar nuclei as a function of 
baryon density rii,. The thick curves lying between the two vertical lines are the results 
in the density region in which the phase with planar nuclei is energetically stable. 

Fig. 10. The critical temperature Tc for the phase with cylindrical nuclei as a function 
of baryon density n^,. The thick curves lying between the two vertical lines are the results 
in the density region in which the phase with cylindrical nuclei is energetically stable. 
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Figure 1: 
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